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We examine the physics capabilities of known flavor neutrino beams from intense muon sources. 
We find that long-baseline neutrino experiments based on such beams can provide precise measure- 
ments of neutrino oscillation mass and mixing parameters. Furthermore, they can test whether the 

fvj ■ dominant atmospheric neutrino oscillations are v^ — > v T and/or v^ — ► v a , determine the z/ M — > v £ 

content of atmospheric neutrino oscillations, and measure v c —* v T appearance. Depending on the 

f^ , oscillation parameters, they may be able to detect Earth matter and CP violation effects and to 

00 ■ determine the ordering of some of the mass eigenstates. 

o 



14.60.Pq, 13.15.+g, 13.35.Bv, 07.77Ka 



I. INTRODUCTION 



a: 

In recent years solar-neutrino and atmospheric-neutrino measurements have provided a growing body of evidence 
for the existence of neutrino oscillations Q . The observed solar neutrino deficit |g,[| can be interpreted as evidence 
for the oscillation of electron neutrinos (v e ) into neutrinos of a different flavor. The recent atmospheric-neutrino 
. __ results from the Super-Kamiokande collaboration B , along with other experiments jq] suggest the oscillation of muon 

S^ ■ neutrinos (v^) into neutrinos of a different flavor, which are dominantly either tau neutrinos (i/,-) or sterile neutrinos 
iy s ). Taken together, these results suggest the mixing of at least three different neutrino types jg-Q]. In addition, the 
LSND collaboration has reported |Hj results from /x + decays at rest that could be interpreted as the first evidence for 
F p — > V e oscillations in an accelerator based experiment. LSND has also reported IllJ results from measurements of 
7r + decays in flight that could be the first evidence for v^ — > v e oscillations. If all of the above reported effects survive, 
then oscillations to a sterile neutrino are required, since there would then be three distinct mass-squared difference 
(Am 2 ) values @. 

The solar-neutrino, atmospheric-neutrino, and LSND results have generated much interest in future accelerator- 
based neutrino oscillation experiments. We can anticipate that an extensive experimental program will be needed 
to firmly establish the existence of neutrino oscillations and to precisely determine all the parameters relevant to 
the phenomenon. The up-coming accelerator-based experiments [|I3] , [l4[ will firmly ground the existence of neutrino 
oscillations, and measure some of the associated parameters. Also, reactor experiments will tightly constrain v e 
disappearance flI5|-|l8|] ; a sensitivity down to ICP 5 cV 2 is expected in Kamland. However, precise measurements of 
neutrino mass-squared differences and the mixing matrix that relates neutrino-flavor eigenstates to neutrino-mass 
eigenstates will almost certainly require a further generation of experiments exploiting higher intensity and/or higher 
quality neutrino beams than currently available |1S|-|21[| . Other goals of such experiments might be fi) to determine 
whether the atmospheric neutrino oscillations are v^ — ► v T , v^ — ► v s , or a mixture of both p2| , p3| , (ii) to detect 
Earth matter effects on neutrino oscillations |24|,g5|], (iii) when matter effects are present, to determine the ordering 
of the neutrino mass eigenstates responsible for the oscillation being measured B5H27J] , and (iv) to detect CP and T 
violation if it exists in the lepton sector Plpq pM ■ 



It has been suggested [[19] that higher intensity and higher quality neutrino beams could be made by exploiting 
the very intense muon sources that are currently being developed as a part of ongoing muon collider |32j feasibility 
studies. The muons would be accelerated up to the desired energy, and injected into a storage ring consisting of two 
long straight sections joined together by two arcs. Muons that decay in the straight sections would form neutrino 
beams consisting of 50% v^ and 50% V e if negative muons are stored, and 50% v e and 50% v ^ if positive muons are 
stored. A compact muon storage ring neutrino source could be tilted downwards at a large angle to enable neutrino 
beams to be directed through the Earth. If the muons from a muon collider muon source are accelerated to energies 
of > 10 GeV and injected into a suitable storage ring, it has been shown that the neutrino fluxes are sufficient to 
detect hundreds of neutrino charged current (CC) interactions per year in a reasonable size detector on the other side 
of the Earth ^. 

The muon storage ring neutrino source idea has led to several recent workshops [ p3| , p4| . A number of pa- 
pers |l9| , po| , p5|j3^ ] have discussed the physics potential of this new type of neutrino facility. In addition, a preliminary 
muon storage ring neutrino source design study has been made J37J, and further studies are planned. The evolution 
of the existing accelerator complexes at Fermilab |38],[39| and CERN Q] towards a muon collider with a muon storage 
ring neutrino source has also been considered. Much further work will be required to develop a realistic design before 
the first muon storage ring neutrino source can be proposed. 

In this paper we consider the physics potential of muon storage ring neutrino sources that are being considered 
in explicit upgrade scenarios for the Fermilab accelerator complex. Two different geometries are considered: (i) a 
neutrino source at Fermilab pointing toward the Soudan mine in Minnesota (L — 732 km), and (ii) a neutrino source 
at Fermilab pointing toward the Gran Sasso underground laboratory in Italy (baseline length L — 7332 km) . 

The paper is organized as follows. In Section II the characteristics of neutrino beams from muon storage ring 
sources are discussed and the basic oscillation formulas are presented. The role of muon storage ring neutrino sources 
in exploring the parameters associated with v^ — > v T and v^ — > v s oscillations is considered in Section III. Section IV 
discusses v e — > v^ oscillations, including matter effects and the possibility of detecting CP violation. In Section V we 
discuss v e — * v T oscillations. Finally, a summary is given in Section VI. 

II. NEUTRINO BEAMS FROM A MUON STORAGE RING 

The neutrino fluxes and interaction rates at a given location downstream of a muon storage ring neutrino source 
depend upon the number of stored muons, the beam divergence within the neutrino-beam forming straight section, 
and the energies and polarization of the decaying muons. In this section we describe the storage ring parameters, 
discuss the muon decay kinematics, present the calculated neutrino fluxes and interaction rates at L = 732 km and 
7332 km, and discuss the matter effects for v e — > is^, v e — + v T , and v^ — * v 8 oscillations. 

A. Storage ring parameters 

A preliminary design study for a muon storage ring neutrino source at Fermilab is described in Ref. |37| . The 
neutrino source consists of: 

(i) An upgraded high-intensity proton source |3q] that cycles at 15 Hz, and delivers 12 bunches per cycle, each 
containing 2.5 x 10 12 protons at 16 GeV. 

(ii) A pion production target, pion collection system, and decay channel. The 16 GeV protons interact in the target 
to produce per incident proton approximately 0.6 charged pions of each sign captured within the decay channel. 
At the end of a 50 m long decay channel on average 0.2 muons of each charge are produced for each proton 
incident on the production target. Hence there would be about 6 x 10 12 muons of the desired charge at the end 
of the decay channel per accelerator cycle, and therefore 9 x 10 20 muons per operational year (10 7 sees). 

(iii) A muon cooling channel that captures the muons exiting the decay channel into bunches with rms lengths 
a z = 1.5 m, and a normalized transverse omittance e^r ~ 0.017 m-rad. The cooling channel reduces the 
omittance to e^r ~ 0.005 m-rad. At the end of the cooling channel there would be about 5.4 x 10 12 muons of 
the desired charge per accelerator cycle, contained within 12 bunches with rms lengths a z = 2 m, a mean muon 
energy of 230 MeV, and an energy spread given by <te/E ~ 0.2. Hence, there would be 8.1 x 10 20 cooled muons 
per operational year. 



(iv) An acceleration system that captures the muons from each long bunch exiting the cooling channel into 16 short 
bunches with a z ~ 0.9 cm, and accelerates the muons to 10 GeV. About 60% of the muons survive the capture 
and acceleration stage. Hence, about 5 x 10 20 muons are accelerated to 10 GeV per operational year. 

(v) A muon storage ring with two long straight sections. The circumference of the ring is 448 m, and the length of 
the neutrino beam-forming straight section is 150 m. Hence about 33% of the muon decays contribute to the 
neutrino beam, and there are 1.6 x 10 20 neutrinos and 1.6 x 10 20 antineutrinos in the beam per operational year. 

Note that in the scenario described above the proton source provides approximately one-third of the beam power 
of the source required for a muon collider. The fluxes and interaction rates discussed in the following sections should 
be multiplied by a factor of 3.3 to obtain the results that would correspond to a muon storage ring neutrino source 
that utilizes the full muon collider muon source. 

B. Fluxes and Interaction Rates 

In the muon rest-frame the distribution of muon antineutrinos (neutrinos) from the decay « — ► e ± + v e (y e ) + 
Vn (ufj.) of polarized muons is given by the expression plj 

£k- f> - 2.) =F(1 -*><»•]. M 

where x = 2E u /m^, 9 is the angle between the neutrino momentum vector and the muon spin direction, and m^ is 
the muon rest mass. The corresponding expression describing the distribution of electron neutrinos (antineutrinos) is 

d 2 N v Ylx 2 . .„ 

-^ = — [ ( l-*)T(l-*)costf]. (2) 

Thus, the neutrino and antineutrino energy- and angular- distributions depend upon the parent muon energy, the 
decay angle, and the direction of the muon spin vector. For an ensemble of muons we must average over the polarization 
of the initial state muons, and the distributions become 

d 2 N 2r 2 

-^ oc — [(3 - 2x) T (1 - 2x)P M cos 9} , (3) 

and 

d 2 N„ Ylx 2 s , 

[{1-x)t(1-x)P ii cos9) , (4) 



dxdVL 4tt 

where P M is the average muon polarization along the chosen quantization axis, which in this case is the beam direction. 
Using Eqs. || and [|, the calculated u^, u e , V^, and V e fluxes at a muon storage ring neutrino source were presented 
in Ref. |19[| as a function of L, P^, and the energy of the stored muons (E^). For unpolarized muons, these results 
can be summarized by the formula 

$=^ = ^ (5) 

dAdt 47ri 2 7 2 (l-/3cosa) 2 ' v ; 

where no is the number of neutrinos (or antineutrinos) per unit time in a given beam, 7 = E^/m^, (3 = p^/E^, 
a is the angle between the beam axis and the direction of interest, and dA is the differential area at the detector. 
In practice the muon beam within the storage ring will have a finite divergence which must be taken into account. 
Current design studies suggest a typical divergence of ~ 1 mr within the beam-forming straight section. The neutrino 
flux at L = 732 km is shown versus a in Fig. [j] for several muon energies. In Fig. g the neutrino flux is shown versus 
the perpendicular distance from the beam axis for several values of L and E^. As Fig. shows, the flux is fairly 
uniform near the axis, and, provided I/7 > 1 mr, begins to fall off only for distances of order L/7 away from the 
beam. In the extreme forward direction cos a ~ 1, and using the approximation /3 ~ 1 — l/(27 2 ), the neutrino flux 
may be approximated as 



These calculated fluxes are summarized in Table | for the baseline lengths L = 732 km and L = 7332 km, and the 
storage ring parameters described Section II. A. 

Equations || and [| can also be used to calculate the distribution dN u /dE v of v^ and v e energies in the beam using 
the approximation that, in the laboratory frame, x — > E v j ' E^ for forward neutrinos from high energy muon decays 
(Efj, >> rri/j,). The calculated v^ and v e spectra are shown in Fig. 0. Note that the v^ spectrum peaks at the stored 
muon beam energy, whereas the v e spectrum peaks at two-thirds of the stored muon beam energy. Also, since $ oc E 2 
and a oc E^, the total event rate is proportional to £?. 

The CC neutrino and antineutrino interaction rates can be calculated using the approximate expressions for the 
cross sections 



a vN ~ 0.67 x 1CT 38 cm 2 x E v (GeV) 



(7) 



and 



a VN ~ 0.34 x 10~ 38 cm 2 x E^ (GeV) 



(8) 



The modification in the linear energy dependence due to the W propagator can be neglected for the neutrino energies 
we consider. Using Eqs. ffl and and assuming no neutrino oscillations, the CC interaction rates corresponding 
to the fluxes in Table B are shown in Table O. Given the distribution of neutrino and antineutrino energies in the 
beam (Fig. ||), and the linear dependence of the CC cross sections with energy (Eqs [?] and g), the predicted energy 
distributions for the resulting CC interactions are shown in Fig. 0. Note that the v^ CC event energies peak at 
the stored muon beam energy, whereas the v e event energies peak at about three-quarters of the stored muon beam 
energy. For v T and v T CC cross sections, we use the ratios acc{v r N) / acciv^N) and <Jcc{vTN)/<Jcc{VfiN) given 
in Ref. J42|, multiplied by the cross sections in Eqs. |7] and |[ If the v T cross sections from Ref. [[43| are used, the 
predicted event rates are 20-30% higher; we use the more conservative values of Ref. J42f . 

Event rates for MINOS are compared to those for a Fermilab to Soudan experiment using a muon storage ring in 
Table III. Since for a muon storage ring neutrino source the total event rate increases as E 3 increasing the muon 
energy dramatically improves the usefulness of such a machine. From this comparison we conclude: (i) for electron 
neutrino beams, even a very low energy muon storage ring might be interesting, and (ii) for muon neutrino beams, 
rates become quite large compared to conventional neutrino beams for E^ > 20 GeV. 

We next consider the distribution of energies of the charged leptons produced in the neutrino and antineutrino CC 
interactions downstream of a muon storage ring neutrino source. It is useful to define the scaling variable 



y 



l-Ei/E v , 



(9) 



where Ei is the energy of a charged lepton produced the CC interaction. The differential neutrino and antineutrino 
cross sections are given approximately by 



da v 
dy 



and 



oc E v 



1 



(i-y) 2 



(10) 



day 

dy 



OC Err 



{i-yf + T 



(ii) 



These formulas neglect the Q 2 dependence of the structure functions. If unpolarized positive muons are stored in 
the muon storage ring neutrino source, the F p CC interaction rate in a distant detector is described by the doubly 
differential distribution obtained from Eqs. [I] and [ll| 



d 2 K, 



cc 



dxdy 



oc x 3 (3 - 2x) 



(l-2/) 2 + 



(12) 



where, for a high energy muon storage ring, x — > E v / Ey, in the laboratory frame. It is convenient to define the 
normalized charged lepton energy 

z = Et/Ey . (13) 

At high energies (1 — y) ~ z/x. Integrating Eq. |l2j over x, the charged lepton energy spectrum is given by 



dz 
yielding 

dN, 



^«jf *«•<»- *)[(:) 



z\ 2 1 
5 



(14) 



oc 1 + 20z 2 - 32z 3 + llz 4 . (15) 



dz 

The corresponding expression for the electron energy spectrum arising from the CC v e interactions is 

dN P - 



dz 



oc 5 + 6z 2 - 32z 3 + 21z 4 . (16) 



The expressions for the \x and e + spectra resulting from CC interactions when unpolarized negative muons are stored 
in the muon storage ring neutrino source are respectively 

—f- oc 5 + 4z 2 - 16z 3 + 7z 4 , (17) 

dz 

and 

^- oc 1 + 30z 2 - 64z 3 + 33z 4 . (18) 

dz 

The charged lepton spectra computed from Eqs. |l5Hl8] are shown in Fig. g. 

C. Interaction Rates with Neutrino Oscillations 

For a given neutrino flavor, neutrino oscillations will modify the neutrino flux at a distant detector, and hence 
the associated charged current interaction rates. Within the framework of two-flavor vacuum oscillations, the flavor 
eigenstates v a and v$ are related to the mass eigenstates Vi and i/j by 

v a = ViCOs9 — VjSW.6 , (19) 

Vj3 = Vismd + Vj cos 6, (20) 

where 9 is the mixing angle. The probability that, while traversing a distance L in vacuum, a neutrino of type a 
oscillates into a neutrino of type /3 is given by 

P(> Q -> vp) = sin 2 (26») sin 2 (1.267Am 2 l i/^) , (21) 

where Atti 2 ^ = m 2 — to 2 is measured in eV 2 /c 4 , L in km, and the neutrino energy E v is in GeV. The neutrino 
oscillation length in vacuum Ly is given by 



Lv = -j-^ ■ (22) 

The first maximum in the oscillation probability occurs when L = Ly/2. The values of Ara| that correspond to this 
oscillation maximum are shown in Fig. |6j as a function of E v for for three baseline lengths. Note that if Am^- is small, 
short baseline lengths require low neutrino energies to probe the oscillation maximum. It is useful to define 

Am 2 , L Am\L , x 

V = —f~ = ~^- x . (23) 

The first maximum in the oscillation probability occurs when r\ = 1.24a;. The modulation of the v e and v^ CC 
interaction spectra for neutrinos originating from a muon storage ring neutrino source is shown in Fig. u\ as a function 
of r\. The distributions of the v a CC interaction energies exhibit peaks whose locations are very sensitive to r\ (and 
hence Am 2 J provided T) ~ 1. 



If there are three neutrino types, then the flavor eigenstates are related to the mass eigenstates by a 3 x 3 unitary 
matrix El 




C12C13 



-iS 



S12C13 



S12C23 - Ci 2 S 2 3Sl3e C12C23 - Sl2S23«13e 



iS 



S12S23 - Cl2C 2 3Sl 3 e 



-C12S23 - Si 2 C23Sl 3 e 




(24) 



where 



cos 6ij and sy = sin 9ij . If the neutrinos are Majorana, there are two extra phases, but these do not affect 



oscillations |4jJ. If Am^ is responsible for solar neutrino oscillations and Am^ for atmospheric neutrino oscillations, 

then lAmlxl < |Am| 

1), appropriate for long-baseline experiments, are 



jil ^ |<-v/;(. 32 |- The resulting vacuum oscillation probabilities for the leading oscillation (1.267|Am 2 1 |_L/E'„ <c 



P(v e -> Vp) = sin 2 023 sin 2 20 13 sin 2 (1.267 Aml 2 L/E u ) , 
P(v e -» v T ) = cos 2 6» 23 sin 2 29 13 sm 2 (1.267 Amj 2 L/E„) , 
P(t/f, -> i/ T ) = cos 4 6»i3 sin 2 26»23 sin 2 (1.267Am§ 2 L/£; j/ ) . 



(25) 
(26) 
(27) 



These expressions are analogous to two-flavor vacuum oscillation probability in Eq. BU except that each oscillation 
channel has a distinct amplitude that depends on the neutrino mixing parameters. 

D. Matter effects 

Electron neutrinos can elastically forward scatter off the electrons in matter via the charged current interaction [B4| . 
When v e oscillates into either v^ or v T , this introduces an additional term in the diagonal element of the neutrino 
flavor evolution matrix corresponding to v e — > v e . It is useful to define the characteristic matter oscillation length La 
as the distance over which the phase of the v e wavefunction changes by 27r |25| : 



u 



2tt 



1.63 x 10 4 



V2G F N e p(g cm" 3 ) Z/A 



km 



(28) 



where N e = pN^Z/A is the electron density in matter of density p, Z/A is the average charge to mass ratio for the 
electrically neutral matter, and No is Avogadro's number. Note that, unlike the vacuum oscillation length Ly, the 
characteristic matter oscillation length Lq is independent of E v . For ordinary rock (p ~ 3 g cm~ 3 and Z/A = 0.5) 
Lq rs 10 4 km. Matter effects are negligible when Lo > L or Lo > Ly. However, matter effects can be appreciable 
for trans-Earth experiments, where L ~ Lo, provided Lq < Ly which is satisfied if Am 2 (eV 2 /c 4 ) < E(GcV)/3000. 
In practice, for a trans-Earth experiment, if E^ is greater than a few tens of GeV the muon storage ring becomes 
too large to have it tilted at a large angle while preserving relatively long straight sections between the arcs. Hence, 
matter effects will only be important for a trans-Earth experiment if Am 2 < 0(1CP 2 ) cV 2 /c 4 . The tilt angle, distance, 
and average electron density for some representative long-baseline experiments from the Fermilab site are given in 
Table |y|. For distances longer than about 10000 km the path goes through part of the core (the core diameter is 
6960 km and the Earth's diameter is 12742 km). 

The effective v e oscillation length in matter L m depends on Lq and Ly. Defining the ratio R m = L m /Ly, it can 
be shown that I25|J47|1 : 



Rn 






Ly 

Lo 



1U 



cos 20 



-1/2 



(29) 



and 



P(v & -> v x ) = R 2 m sin 2 (29) sin 2 



1.267 Am]; L 



Rn 



E„ 



(30) 



For antineutrinos, N e in Eq. ^8| changes sign; hence matter effects are different for neutrinos and antineutrinos if 
cos 29 7^ 0. If cos29Am 2 i > (i.e., v e is more closely associated with i/, than Vj), oscillations of v e are enhanced and 
oscillations of v e are suppressed when Ly ~ Lq; if cos26*Am 2 j < 0, the situation is reversed. Therefore, matter effects 
may distinguish between the two different mass orderings when cos 29 =/= 0. 

If Ly/Lo = cos 29, the oscillations in matter have maximal mixing [ p5|j . Figure pi shows the values of Am 2 /E u versus 
sin 2 29 that give maximal mixing in the mantle and core of the Earth; similarly, Fig. g shows values of E v versus sin 2 29 



that give maximal mixing in the Earth for Am 2 = 3.5 x 10~ 3 eV 2 /c 4 , the value favored by the Super-Kamiokande 
atmospheric 708-day data [Q . 

Note that both the oscillation amplitude and the oscillation length depend on R mi and hence on the density p. 
Fortunately, outside of the core the density profile of the Earth is well known, and is described by the Preliminary 
Reference Earth Model p9| . Density profiles along a selection of chords passing through the Earth are shown in 
Fig. [Kj. For L = 732 km most of the pathlength is in rock with p ~ 2.5 g cm~ 3 . For L = 7332 km most of the 
pathlength is in higher density matter with p > 4 g cm -3 . In some early studies [E5LpO| matter effects were computed 
using the average density p av along the traversed path. This gives a reasonable estimate in many cases, but can 
introduce significant errors in the calculation of P(v e — > v x ) in some regions of parameter space. To illustrate this, 
Fig. [Li] shows P(y e — > v x ), for x = p or r, versus p calculated using Eq. [3(j for 10 GeV neutrinos propagating 7332 km. 
The calculation assumes sin 2 29 = 1, and the results are shown for three different values of Am 2 . For L = 7332 km 
the average density is p av — 4.2 g/cm 3 . Consider the Am 2 ~ 0.001 eV 2 /c 4 curve, for which Ly = 24800 km. If p av 
is used we obtain L = 7810 km and R m = 0.30. The resulting oscillation probability is P(p av ) = 0.002. However, 
nearly all of the pathlength is through matter with densities that correspond to oscillation probabilities significantly 
higher than P(p av )- If the true density profile is used instead of p av we obtain P = 0.014. Thus, in this example, 
using p av leads to a result which underestimates the oscillation probability by a factor of 7. Therefore, in the following 
sections we do not use p av to calculate P(y e — > v x ), but integrate using the density profile given by the preliminary 
Earth model. Extensive analyses have been made of oscillation effects involving transmission through the Earth's 
mantle and core |27|,|l|-||] . 

For v^ oscillations to v Tl there are no matter effects for simple two-flavor oscillations (although there may be small 
matter effects on v^ — > v T for three flavors [p5[). For v^ — > v s oscillations, however, muon neutrinos elastically forward 
scatter off the quarks and electrons in matter via the neutral current interaction, whereas sterile neutrinos do not. 
The matter oscillation length in this case can be found by replacing N e in Eq. |28| by ~N n /2 (by N n /2 for v^ — v v s 
oscillations) [ p6| ; Z in Eq. pq is then replaced by Z — A (A — Z). If the mixing of atmospheric v^ is nearly maximal 
(as suggested by atmospheric neutrino data), then Eqs. E3 and au imply that the amplitude of v^ — ► v s and v^ — ► v s 
oscillations will both be suppressed by matter effects if Ly ^ Lo, and by the same amount since cos 28 ~ 0. 

III. u„ -» v T AND i/„ -> us OSCILLATIONS 

The currently favored explanation for the Super-Kamiokande atmospheric neutrino results is that muon neutrinos 
are oscillating primarily into either tau neutrinos or sterile neutrinos, with the oscillation parameters given by sin 29 ~ 
1 and Am 2 in the approximate range 0.002 eV 2 /c 4 to 0.006 eV 2 /c 4 P,H- Searches for neutrino oscillations at 
accelerators are based on either the disappearance of the initial neutrino flavor, or the appearance of a neutrino flavor 
not originally within the neutrino beam. If the current interpretation of the Super-Kamiokande data is correct, the 



next generation of approved long baseline accelerator experiments [TL3 14 should confirm the existence of neutrino 
oscillations and make the first laboratory measurements of the oscillation parameters. This will happen in the period 
before a first muon storage ring neutrino source could be built. It is likely therefore that the main atmospheric neutrino 
oscillation physics goals in the muon storage ring neutrino source era would be to make very precise measurements 
of the oscillation parameters, and determine whether there is a small v^ — » v s (y^ — > v T ) component in a dominantly 
Vp — ► v T (v^ — > v s ) signal. CC measurements can distinguish between v^ — > v T and v^ — ► v s oscillations in two ways: 
(i) direct appearance of taus for v^ — > v T |57J and (ii) a different v^ — » v s transition probability due to matter effects 
when compared to v^ — > v r . NC/CC measurements of ir° production can also distinguish v^ — » v s from v^ — » v r |p8| . 

A. Fermilab — > Gran Sasso 

Consider first a 10 GeV muon storage ring at Fermilab with the neutrino beam-forming straight section pointing 
at Gran Sasso (L = 7332 km). The predicted CC event rates are listed as a function of E^ and Am 2 in Table [V|. The 
predicted v T appearance event rates arising from v^ — * v T oscillations with sin 28 = 1 and Am 2 in the range from 
0.002 eV 2 /c 4 to 0.006 eV 2 /c 4 are only of order 1 event per kt-yr. This event rate is too low to enable a precise v T 
appearance measurement in, for example, a 1 kt hybrid emulsion detector. However, an oscillation measurement could 
be made by a f M disappearance experiment. In the absence of oscillations about 220 v^ CC events would be expected 
per year in a 10 kt detector, with an event energy distribution peaking at ~ 10 GeV. In the presence of oscillations 
the predicted v^ CC event rate (Table [v]) and the corresponding event energy distribution (Fig. 12) are both sensitive 



to Am 2 (provided r\ ~ 1). However in general, for a given Am 2 , the average oscillation probability will depend on 
whether v^ — > v T or v^ — > v a or a mixture of the two oscillations is being observed. Hence, to avoid a large uncertainty 



in the extracted value of Am 2 due to an uncertainty in which oscillation mode is being observed, it is desirable that 
a trans-Earth baseline experiment is sensitive to both v^ disappearance and v T appearance. We conclude that the 
muon storage ring beam energy needs to be higher than 10 GeV to facilitate a v T appearance measurement. Consider 
next a 20 GeV muon storage ring at Fermilab with the neutrino beam-forming straight section pointing at Gran Sasso 
(L = 7332 km). The v r CC rate corresponding to the favored region of v^ — > v T parameter space is now sufficiently 
large to facilitate a v T appearance measurement. In the absence of oscillations about 1900 v^ CC events would be 
expected per year in a 10 kt detector. In the presence of v^ — ► v T oscillations with sin 2 29 = 1 and Am 2 = 0.002 eV 2 /c 4 
(0.006 eV 2 /c 4 ) only 370 (1300) v^ CC events per 10 kt-yr would be expected, together with 43 (14) v T CC events 
per year in a 1 kt detector. The corresponding event rates for v^ — ► v s oscillations are 1400 (1200) v^ CC events 
per 10 kt-yr and no v T CC events. With these trans-Earth event rates at a 20 GeV storage ring a v^ disappearance 
measurement would enable the oscillation probability to be measured with a precision of a few percent and, for a 
known mixture of v^ — > v r and v^ — > v s oscillations, would enable Am 2 to be determined with a statistical precision 
of a few percent. In addition, a v T appearance measurement would enable the z/„ — * v T oscillation probability to 
be measured with a statistical precision which depends on Am 2 , but is typically about 20%. Finally, a comparison 
of the appearance and disappearance results would enable a v^ — > v a contribution at the few times 10% level to be 
observed in a predominantly v^ — > v r signal, or an approximately 10% v^ — ► v r contribution to be measured in a 
predominantly v^ — > v s signal. 

B. Fermilab —* Soudan 

Consider a 10 GeV muon storage ring at Fermilab with the neutrino beam-forming straight section pointing at 
Soudan (L — 732 km). Predicted event rates are listed in Table [V| In the absence of oscillations about 22200 v^ CC 
events would be expected per year in a 10 kt detector. In the presence of either v^ — > v T or v^ — ► v s oscillations with 
sin 2 20 = 1 and Am 2 = 0.002 eV 2 /c 4 (0.006 eV 2 /c 4 ) only 20500 (11700) v^ CC events per 10 kt-yr would be expected. 
In addition, for v^ — > v T oscillations 20 (150) v T CC events per year would be expected in a 1 kt detector. Note that the 
shape of the v T CC interaction energy distribution (Fig. [13|) exhibits some dependence on Am 2 , although a detector 
with good energy resolution would be needed to exploit this dependence. We conclude that a storage ring with a 
beam energy as low as 10 GeV might be of interest for a Fermilab — ► Soudan experiment, particularly if Am 2 is at 
the upper end of the currently favored region. In this case the average oscillation probability could be measured with 
a statistical precision of better than 1% from the disappearance measurement, and the resulting Am 2 be determined 
with a precision of ~1%. Note that, since matter effects are small, there is no additional uncertainty on Am 2 arising 
from an imprecise knowledge of the oscillation mode {v^ — > v T or v^ — ► v s ). Finally, with Am 2 = 0.006 eV 2 /c 4 , a 
comparison of the v^ disappearance and v T appearance measurements would enable a few percent contribution from 
v^ — ► v T oscillations to be measured in a predominantly v^ — > v s signal, or a few times 10% v^ — > v s contribution to 
be measured in a predominantly v^ — v v T signal. 

In Fig. H4 single v T CC event per kt-yr contours in the (sin 2 20, A?n 2 )-plane are shown for L = 732 km as a 
function of E^. The sensitivity of the v T appearance measurement increases with increasing muon storage ring beam 
energy. Consider the sensitivity for a 50 GeV muon storage ring. The dependence of the v T CC rates on sin 2 29 and 
Am 2 is shown in Fig. |l5j. The predicted event rates range from several hundred events per kt-yr at the lower end 
of the Super-Kamiokande allowed region of parameter space, to several thousand events per kt-yr at the upper end 
of the favored region (Table |V|). Hence the v^ — ► v T oscillation probability could be measured with a precision of 
a few percent, corresponding to a determination of Am 2 with a precision of a few percent. However, the predicted 
oscillation probabilities are small, and the statistical precision of the disappearance measurement is less than the 
corresponding precision for the E^ = 10 GeV case. Note that about 2 million v^ CC events would be expected per 
year in a 10 kt detector. The sensitivity of the disappearance measurement may therefore be limited by systematics. 
This needs further study. We conclude that, although the statistical sensistivity of the v T appearance measurement 
improves with increasing E^, the optimum choice of muon beam energy depends upon whether Am 2 is in the lower 
or upper part of the favored region of parameter space, and on at what level the higher energy experiments would be 
limited by systematics. 

C. L = 7332 km versus 732 km 

The relative performance of experiments at long [L = 732 km) and very long [L = 7332 km) baselines depends on 
both Am 2 and the oscillation mode (y^ — > ^ T or v^ — > v s ). Based on the results summarized in Table |V| we find that: 



Case 1: 0.002 < Am 2 < 0.004 eV 2 /c 4 and v^ — > v T oscillations dominate. An L — 732 km experiment with 
Ept = 10-20 GeV, or an L = 7332 km experiment with E^ = 20 GeV both seem interesting. A higher energy 
L = 732 km experiment (e.g. E^ = 50 GeV) would provide a higher statistics measurement of the v^ — > j^ r 
oscillation probability, but the oscillation probability is small and the measurement precision may be dominated 
by systematics. 

Case 2: 0.004 < Am 2 < 0.006 eV 2 /c 4 and v^ — > ^ r oscillations dominate. An i = 732 km experiment with B^ = 10- 
20 GeV would provide a good v^ disapearance measurement and a good v T appearance measurement. With a 
larger E^ (e.g. 50 GeV) the unoscillated neutrino CC rate becomes large and the oscillation measurements may 
be limited by systematics. The trans-Earth baseline option with E^ = 20 GeV yields less statistical precision 
than the L = 732 km option with E^ = 20 GeV. 

Case 3: v^ — > v s oscillations dominate. In this case a very sensitive v^ disapearance measurement is important. A 
low energy (10 GeV) storage ring pointing at Soudan seems attractive, particularly if Am 2 is in the upper half 
of the favored region. An L = 7332 km baseline experiment seems to be less attractive than an L = 732 km 
baseline experiment for large Am 2 , but is complementary to the shorter baseline experiment for small Am 2 , 
and would enable matter effects to be measured. 

IV. u e -> ^ M OSCILLATIONS 

If confirmed, the LSND neutrino oscillation results necessarily imply the existence of v e — > v^ oscillations, with 
0.3 < Am 2 LSND < 2.0 eV 2 /c 4 . Oscillation limits from the Karmen experiment |39[ are no longer in conflict with the 
LSND effect. In addition, the solar neutrino results can also be interpreted in terms of v e — > v^ and/or v e — > v T 
oscillations either within the framework of the MSW effect with Am 2 un ~5x 10~ 5 eV 2 /c 4 and sin 2 29 ~ 0.8 (large 
angle solution) or Am 2 un ~5x 10 -6 eV 2 /c 4 and sin 2 29 ~ 5 x 10~ 3 (small angle solution) p0-p2|, or within the 
framework of vacuum oscillations |33[ with Am 2 un ~ 0.65 — 8.6 x 10 -10 eV 2 /c 4 [|[^0|. In the context of three neutrinos, 
although the primary oscillation of atmospheric neutrinos is f„ — > v T , a small v„ —* v e component is not excluded. For 
Am 2 tm ~ 3.5 x 10 -3 eV 2 /c 4 , a measurement of v e — > v^ will place strong limits on sin 2 2 3 sin 2 2^i 3 (see Eq. |25|) fl |[ ■ 

The next generation of approved neutrino experiments include MiniBooNE |54j at Fermilab, and KamLAND [gj at 
the Kamioka Mine, Japan. MiniBooNE, which will consist of a 445 ton fiducial volume liquid scintillator detector 500 m 
downstream of a v^ source, with E v ~ 0.75 GeV, should confirm or refute the neutrino oscillation interpretation of the 
LSND results. KamLAND, which will consist of a 1 kt liquid scintillator detector 140 km and 200 km from reactors in 
Japan, should observe a v e disappearance signal if the solar neutrino results are due to oscillations corresponding to the 
large angle MSW solution. Both MiniBooNE and KamLAND nominally begin operation in ~ 2001. In the following 
we consider the v e — ► v^ physics potential of a muon storage ring neutrino source in the post-MiniBooNE/KamLAND 
era. 

A. Fermilab — » Soudan 

Figure |l6| shows single event per year contours in the (sin 2 29, Am 2 )-plane for v^ CC interactions in a 10 kt detector 
732 km downstream of the muon storage ring neutrino source described in Section II. A, where the v^ arise from v e — > v^ 
oscillations. The contours are shown for several values of E^ from 10 GeV up to 250 GeV. The upper part of the MSW 
large mixing angle solution is within the single-event per year boundary for E^ > 10 GeV. The sensitivity improves 
with increasing E^ . The event rate is shown in the sin 2 29- Am 2 plane in Fig. O for E^ — 250 GeV. To completely 
cover the parameter space corresponding to the MSW large mixing angle solution would require larger muon beam 
currents than those of the scenario described in Section II. A. It should be noted that the full muon collider front-end 
described in Ref. |32j| , would produce a factor of 3.3 more muons. The resulting neutrino flux would just about enable 
the entire MSW large mixing angle solution to be covered by a 50 GeV storage ring and a few years of running with, 
for example, a 20 kt detector having a detection efficiency of 50%, provided background rates are negligible. The 
small angle MSW solution is out of reach of a muon storage ring neutrino source, as currently envisioned. Hence if 
KamLAND rules out the large angle solution, then solar neutrino oscillations would provide little motivation for a 
Fermilab — + Soudan v e — ► v^ appearance experiment. On the other hand, should KamLAND obtain a positive signal 
for solar neutrino oscillations corresponding to the large angle MSW solution, a Fermilab — > Soudan muon storage 
ring neutrino source experiment to search for v e — > v^ appearance might be feasible, but would clearly be challenging 
and would probably require the full muon collider muon source for several years. If this becomes the only way to 



discover whether the solar neutrino oscillations are due to v e — > v^ or due to v e — > v T or v s , then this experiment 
would certainly be worthy of serious consideration. 

Independent of the KamLAND results, a positive MiniBooNE result would motivate precise measurements of 
v e — ► Vfj, oscillations to test whether or not P(y e — » t^) = P(^ — > v e ) (T invariance), test the oscillation phenomenon 
with a very different L and E u , search for a v e — > ^ s contribution, etc. A Fermilab — ► Soudan j/ e — > t^ experiment 
would enable many thousands of v^ appearance CC interactions (due to LSND-scale oscillations) to be measured (see 
Figs. [H] and [jj]) . Such measurements could also test the existence of v e — + v^ oscillations at the atmospheric Am 2 scale 
(~ 3.5 x 10~ 3 cV 2 /c 4 ), which allows a precise determination of the three- neutrino mixing parameter cos 2 6*23 sin 2 29±s. 

B. Matter effects 

Figure O shows single event per year contours in the (sin 2 29, Am 2 )-plane for v^ CC interactions in a 10 kt detector 
7332 km downstream of the muon storage ring neutrino source described in Section II. A, where the ia. arise from 
v e — ► Vy, oscillations with Am 2 > 0. The contours are shown for E^ = 10, 20, 50, and 250 GeV. Figure |l^ shows the 
corresponding contours for v e — ► P^ oscillations. Compared with the equivalent contours for L = 732 km (Fig. (l6| ) the 
lower event statistics suppress the sensitivity at large Am 2 and small sin 2 29, and matter effects (Eq. £30|) suppress the 
sensitivity at small Am 2 and large sin 2 29. It is interesting to consider whether the matter effects might be measured 
by a trans-Earth muon storage ring neutrino source v e — ► Vy. experiment searching for wrong-sign muons. 

Note that if MiniBooNE confirms the LSND evidence for oscillations, then Am 2 will be too large to result in a 
significant modification of the oscillation probability due to matter effects. However, if KamLAND obtains a positive 
v e disappearance result corresponding to the large mixing angle MSW solution for solar neutrino oscillations, then the 
Am 2 will be in the right range to produce a significant modification of the oscillation probability due to matter. As 
an example, consider a storage ring with E^ = 20 GeV, and assume that sin 2 29 = 1. Then if Am 2 = 1 x 10~ 4 eV 2 /c 4 , 
in the absence of matter effects ~ 11 wrong sign muons (from v e — ► v^ oscillations) would be expectected per 10 kt-yr. 
With matter effects, the number of wrong sign muons is reduced to ~ 0.6 per 10 kt-yr. Hence, an experiment to 
measure matter effects, by comparing the wrong sign muon rates in a long baseline (e.g L = 732 km) experiment with 
the corresponding rates in a trans-Earth (e.g L — 7332 km) experiment might be feasible if the solar neutrino results 
are due to v e — > v^ oscillations with sin 2 29 ~ 1 and Am 2 at the upper end of the region of allowed parameter space 
for the MSW large mixing angle solution. This experiment would clearly be challenging, and would probably need a 
more intense neutrino source than described in Section II. A (for example, the full muon collider muon source utilized 
for several years with a 20 kt detector). 

If there is a small v^ — ► v e and P M — > P e component to atmospheric neutrino oscillations, then v e — > v^ and 
v e — y v u oscillations are also expected; a trans-Earth experiment could test for different matter effects acting on these 
two oscillation channels. Table W\ shows, for a Fermilab — > Gran Sasso experiment (L = 7332 km), the v e — > v^ 
and v e — > v^ oscillation probabilities, the predicted numbers of v^ and P M CC appearance events, and the average 
probability for v e — ► v^ assuming no matter effects. The predictions are tabulated for sin 29 = 0.1 with various E^ 
and values of Am 2 > suggested by the atmospheric neutrino data. If Am 2 < (i.e., the lower mass eigenstate 
is more closely associated with v^ and not v e ), then the oscillation probabilities for neutrinos and antincutrinos in 
Table ylt are reversed. Thus the sign of Am 2 may be determined from these measurements if matter effects are 
present~ |25|J27| l . 

In this paper we have assumed that the neutrino mixing conserves CP; the apparent violation of CP in the v e — ► v^ 
and v e — » Vy, channels would be due to matter effects. Measuring v^ — > v e and P M — ► v e oscillations as well (through an 
electron appearance experiment) would provide important cross-checks and help distinguish between matter effects 
and CP violation. For example, CPT conservation for vacuum oscillations implies P(v e — * Vp.) = P{v^ — > v e ), 
whether CP is violated or not, so that a difference between these two oscillation probabilities would be unambiguous 
evidence for matter effects. Similarly, matter effects are the same for v e — ► v^ and v^ — > v e oscillations, so that 
a difference between these oscillation channels is due to T violation, which in turn implies CP violation under the 
plausible assumption that CPT is conserved. 

V. v e -*■ v T OSCILLATIONS 

If the atmospheric neutrino results are due to v^ — > v T oscillations with a Am 2 ~ 3.5 x 10 -3 eV 2 /c 4 , and the solar 
neutrino results are due to v e — > v^ and/or v e — > v T oscillations with a Am 2 < O(10 -5 ) eV 2 /c 4 , it is tempting to 
conclude that v e — > v T oscillations should exist with an associated Am 2 = O(10 -3 ) eV 2 /c 4 , but with unknown effective 
sin 2 26. For Am 2 > 2 x 10~ 3 cV 2 /c 4 , v e — * v T oscillations with sin 2 29 > 0.2 are already excluded by the CHOOZ 
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reactor v e disappearance experiment, although there is no limit for Am 2 < 10~ 3 eV 2 /c 4 . In a three-neutrino scenario, 
measurements of atmospheric neutrinos indicate sin 2#23 > 0.8 and exclude sin 29i3 > 0.33 for all atmospheric 
Am 2 jjj; v e — > v T oscillations may be detectable if #13 is not too small. Hence, the region of parameter space that 
would be interesting to cover in future v e — > v r oscillation searches is given by Am 2 ~ 3.5 x 10 -3 eV 2 /c 4 and effective 
sin 2 29 < 0.3. 



Consider first a v e — ► v r appearance experiment. Event rates are summarized in Table VI] for the neutrino source 
described in Section II. A, and sin 2 29 — 0.1. For a baseline length of 732 km matter effects are small, and hence the 
oscillation rates for sin 2 29 < 0.1 can be obtained by multiplying the rates in the table with the value of sin 2 20/0.1. For 
a baseline length of 7332 km matter effects are large and enhance the v T appearance rate for Am 2 ~5x 10 -3 eV 2 /c 4 , 
although the numbers of events are not large. On the other hand, the corresponding rates for a 1 kt detector 732 km 
downstream yield more v T appearance events in the region of interest, and for En ~ 50 GeV could yield up to a few 
events per year with sin 2 29 as small as 10 -2 (see Fig. g(j). Single-event per kt-yr sensitivity contours in the (sin 2 29, 
Am 2 )-plane are shown versus the stored muon energy in Fig. |2l]. Decreasing the storage ring energy from 50 GeV to 
20 GeV would increase the minimum sin 2 29 that could be probed by a factor of a few. 

In a three-neutrino scenario, a 1 kt detector 732 km downstream could probe cos 2 #23 sin 2 26>i3 to very low values; 
when combined with v e — > v^ measurements (see Sec. IV) and with future atmospheric neutrino results, a precise 
determination of #23, #13 and Am 2 should be possible. 

VI. SUMMARY 

In this paper we have made a comprehensive study of possible neutrino oscillation measurements in long-baseline 
experiments with an intense neutrino beam from a muon storage ring. Our results can be summarized as follows: 

(i) A neutrino beam from a muon storage ring provides much larger electron neutrino interaction rates than conven- 
tional neutrino beams, and much larger muon neutrino interaction rates if E^ > 20 GeV. 

(ii) If the Am 2 associated with the atmospheric neutrino results is in the lower part of the allowed Super-Kamiokandc 
parameter space (Am 2 < 4 x 10~ 3 eV 2 /c 4 ), a Fermilab — > Gran Sasso trans-Earth experiment with an E^ = 
20 GeV muon storage ring neutrino source would enable v^ disappearance to be established in an accelerator- 
based experiment, and Ato 2 to be measured with a precision of a few percent. Matter effects may also allow 
the differentiation of v^ — > v r and v^ — ► v s oscillations. 

(iii) If the Am 2 associated with the atmospheric neutrino results is in the upper part of the allowed Super- 
Kamiokande/Kamiokande parameter space (Am 2 > 4 x 10~ 3 eV 2 /c 4 ), a Fermilab — > Soudan (or CERN — > 
Gran Sasso) v^ disappearance experiment would allow Am 2 to be measured with a precision of a few percent; 
a v T appearance experiment would enable a search for a small contribution from v^ — > v s within a dominantly 
v^ — ► v T signal, or conversely, a small contribution from v^ — ► v T within a dominantly v^ — > v s signal. 

(iv) If the LSND v^ — > v e oscillation results are confirmed by MiniBooNE, a Fermilab — > Soudan v e — ► v^ experi- 
ment would enable many thousands of wrong-sign muons to be measured, facilitate a precise confirmation and 
measurement of the oscillation parameters, and enable Piy^ — ► v e ) to be compared with P(y e — > v^) for a test 
of T invariance. 

(v) If KamLAND observes v e disappearance with a rate corresponding to the large angle MSW solar neutrino solution 
then a measurement of v e — > v^ appearance might be concievable for a Fermilab — > Soudan experiment if the 
full muon source needed for a muon collider were available for several years of running with a muon storage 
ring neutrino source. It might also be possible to measure matter effects as the electron neutrinos traverse the 
Earth. These measurements are however challenging and need further consideration to assess their viability. 

(vi) If there is a small but non-negligible component of v^ «-> v e oscillations in atmospheric neutrinos, a Fermilab — > 
Gran Sasso measurement of v e — ► v^ and v e — > v^ could perhaps detect matter effects in the oscillation, which 
show up as apparent CP violation. Matter effects and true CP violation can be distinguished by also measuring 
v^ — > v e and v^ — ► v e oscillations: a nonzero P{v e — > v^) — Piy^ — > i> e ) signals T (and hence CP) violation, a 
nonzero P{y e — > v^) — Piy^ — > v e ) signals matter effects, while a nonzero P{v e — > v^) — P{v e — > P M ) can occur 
due to either matter or CP violation. If oscillations of neutrinos are enhanced and antineutrinos suppressed, 
for the reverse situation, 
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(vii) A v e — ► v T appearance search sensitive to Am 2 O(10~ 3 ) eV 2 /c 4 and sin 2 2d down to about 10 -2 seems possible 
with a Fcrmilab — ► Soudan experiment using a > 20 GeV muon storage ring. Higher energies are more desirable 
for v T appearance experiments. 

(viii) In a three-neutrino scenario with different Am 2 scales for the solar and atmospheric oscillations, measurements 
of u e — ► Up and v e — > v T should provide an accurate determination of the two mixing angles associated with the 
leading oscillation. When combined with solar neutrino measurements, the three-neutrino mixing parameters 
would be completely determined, except for a possible CP violating phase. 
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TABLE I. Neutrino and antineutrino fluxes calculated for the baseline lengths L — 7332 km (FNAL — ► Gran Sasso) and 
L — 732 km (FNAL — ► Soudan). In the calculation the neutrino fluxes have been averaged over a circular area with radius 
1 km (0.1 km) for L — 7332 km (732 km). This averaging provides an approximate model for the 0(10%) reduction in the 
neutrino fluxes at the far site due to the finite divergence of the parent muon beam in the straight section of the storage ring. 



parent 


£ M (GcV) 


10 


20 


50 


250 




FNAL -> Soudan 










^ 


v„ (10 12 m,- 2 yv- 1 ) 


0.79 


3.3 


21 


460 


^ 


VI (10 12 m~ 2 yr" 1 ) 


0.79 


3.3 


21 


470 


M+ 


v e (10 12 m- 2 yr" 1 ) 


0.85 


3.3 


21 


470 


M+ 


TV (10 12 mT 2 yr" 1 ) 


0.84 


3.3 


21 


470 




FNAL -s- Gran Sasso 










^ 


v„ (10 10 m- 2 yr" 1 ) 


0.79 


3.3 


21 


470 


^ 


—, (10 10 m- 2 yr" 1 ) 


0.79 


3.3 


21 


470 


M+ 


v„ (10 10 m-' 2 yr" 1 ) 


0.83 


3.3 


21 


470 


M+ 


TV (10 10 m- 2 yr" 1 ) 


0.85 


3.3 


21 


470 



TABLE II. Neutrino and antineutrino CC interaction rates in the absence of oscillations, calculated for the baseline lengths 
L = 7332 km (FNAL — ► Gran Sasso) and L — 732 km (FNAL — > Soudan). In the calculation the neutrino fluxes have been 
averaged over a circular area with radius 1 km (0.1 km) for L = 7332 km (732 km). This averaging provides an approximate 
model for the O(10%) reduction in the neutrino fluxes at the far site due to the finite divergence of the parent muon beam in 
the straight section of the storage ring. 



parent 


E„ (GeV) 




10 


20 


50 


250 




FNAL -> Soudan 












^ 


Vp (per kt-yr) 


2.2 


x 10 3 


1.9 x 10 4 


2.9 x 10 5 


3.1 x 10 7 


M~ 


U2 (per kt-yr) 


9.6 


x 10 2 


8.0 x 10 3 


1.3 x 10 5 


1.4 x 10 7 


M+ 


v e (per kt-yr) 


1.9 


x 10 3 


1.6 x 10 4 


2.4 x 10 5 


2.7 x 10 7 


M+ 


7v (per kt-yr) 


1.2 


x 10 3 


9.5 x 10 3 


1.5 x 10 5 


1.6 x 10 7 




FNAL -> Gran Sasso 












^ 


Vp (per kt-yr) 




22 


1.9 x 10 2 


2.9 x 10 3 


3.1 x 10 5 


M~ 


vl (per kt-yr) 




9.5 


80 


1.3 x 10 3 


1.4 x 10 5 


^ + 


v e (per kt-yr) 




20 


1.6 x 10 2 


2.4 x 10 3 


2.7 x 10 5 


M + 


7v (per kt-yr) 




12 


94 


1.5 x 10 3 


1.6 x 10 5 



TABLE III. Muon neutrino and electron antineutrino CC interaction rates in the absence of oscillations, calculated for 
baseline length L — 732 km (FNAL — ♦ Soudan), for MINOS using the wide band beam and a muon storage ring with 
Ef, = 10, 20, 50 and 250 GeV. 



Experiment 




(GeV) 


(GeV) 


N(*/ M CC) 
(per kt-yr) 


N(^ e CC) 
(per kt-yr) 


MINOS 


Beam 






Low energy 

Medium energy 

High energy 


3 
6 
12 


- 


458 
1439 
3207 


1.3 
0.9 

0.9 


Muon ring 


E» (GeV) 






10 
20 
50 
250 


7 
14 
35 
175 


6 
12 
30 
150 


2200 

19000 

2.9 xlO 5 

3.1 xlO 7 


960 

8000 

1.3 xlO 5 

1.4xl0 7 
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TABLE IV. Tilt angle, distance, and average electron density for long-baseline experiments with neutrino source at Fermilab. 



Detector 
site 


Tilt angle 
(degrees) 


L 

(km) 


(No/cm*) 


J-'core 

(km) 


V*e/ mantle 

(N Q /cm 3 ) 


\-f » e) core 

(AV cm 3 ) 


Soudan 


3.3 


732 


1.67 





1.67 


- 


Gran Sasso 


35 


7332 


2.09 





2.09 


- 


Kamioka 


46 


9160 


2.26 





2.26 


- 


South Pole 


67 


11700 


3.39 


4640 


2.33 


4.99 



1G 



TABLE V. Muon neutrino CC interaction rates calculated for a distance L downstream of a storage ring containing 



unpolarized negative muons. The event rates assume v^ — > u T oscillati 


ons with sin 28 


= 1. Results for u 


_i — > u s oscillations 


are shown 


in parentheses 


for L — 7332 km; for L — 732 km matter effects are small and the results for u^ 


— ► u T and Uy, — ► u s 


are very similar. Also shown are the tau neutrino CC interaction rates due to u^ — > u T 


Dscillations. 




E» 


L 


Am 2 


<^> 


N(k m CC) 


Significance 


N(u T CC) 


(GeV) 


(km) 


(eV 2 /c 4 ) 


J-V - * v t { V V —* V s) 


(per 10 kt-yr) 




(per kt-yr) 


10 


7332 





(0) 


220 (220) 


0(0) 





10 


7332 


0.002 


0.53 (0.18) 


100 (180) 


12cr (3.0<t) 


1.9 


10 


7332 


0.003 


0.29 (0.36) 


160 (140) 


5.2a (6.7a) 


0.7 


10 


7332 


0.004 


0.63 (0.64) 


80 ( 80) 


15a (16a) 


2.1 


10 


7332 


0.005 


0.56 (0.48) 


100 (120) 


12a (9.9a) 


1.9 


10 


7332 


0.006 


0.38 (0.37) 


140 (140) 


7.1a(7.0a) 


1.2 


20 


7332 





(0) 


1900 (1900) 


0(0) 





20 


7332 


0.002 


0.80 (0.27) 


370 (1370) 


78a (14a) 


43 


20 


7332 


0.003 


0.81 (0.25) 


360 (1400) 


80a (12a) 


46 


20 


7332 


0.004 


0.53 (0.18) 


880 (1530) 


34a (8.6a) 


30 


20 


7332 


0.005 


0.30 (0.21) 


1310 (1480) 


15a (10a) 


15 


20 


7332 


0.006 


0.29 (0.36) 


1320 (1190) 


15a (20a) 


14 


50 


7332 





(0) 


28800 (28800) 


0(0) 





50 


7332 


0.002 


0.26 (0.10) 


21300 (25900) 


52a (18a) 


350 


50 


7332 


0.003 


0.48 (0.18) 


15000 (23700) 


110a (33a) 


660 


50 


7332 


0.004 


0.67 (0.24) 


9500 (22000) 


200a (46a) 


950 


50 


7332 


0.005 


0.80 (0.27) 


5700 (21000) 


310a (54a) 


1160 


50 


7332 


0.006 


0.86 (0.28) 


4200 (20800) 


380a (55a) 


1260 


10 


732 








22200 








10 


732 


0.002 


0.08 


20500 


12a 


20 


10 


732 


0.003 


0.16 


18700 


26a 


43 


10 


732 


0.004 


0.26 


16500 


45a 


74 


10 


732 


0.005 


0.37 


14100 


69a 


110 


10 


732 


0.006 


0.47 


11700 


98a 


150 


20 


732 








185000 








20 


732 


0.002 


0.020 


181000 


8.6a 


89 


20 


732 


0.003 


0.044 


177000 


19a 


200 


20 


732 


0.004 


0.076 


171000 


34a 


350 


20 


732 


0.005 


0.11 


164000 


52a 


530 


20 


732 


0.006 


0.16 


156000 


74a 


750 


50 


732 








2.88 x 10 B 








50 


732 


0.002 


3.3 x 10" 3 


2.88 x 10 6 


5.6a 


410 


50 


732 


0.003 


7.4 x 10 -3 


2.86 x 10 6 


13a 


910 


50 


732 


0.004 


0.013 


2.84 x 10 6 


22a 


1600 


50 


732 


0.005 


0.020 


2.83 x 10 6 


35a 


2500 


50 


732 


0.006 


0.029 


2.81 x 10 6 


50a 


3600 


250 


732 








3.11 x 10 s 








250 


732 


0.002 


1.4 x 10~ 4 


3.11 x 10 8 


2.5a 


3000 


250 


732 


0.003 


3.2 x 10~ 4 


3.11 x 10 8 


5.7a 


6800 


250 


732 


0.004 


5.8 x 10~ 4 


3.11 x 10 8 


10a 


12100 


250 


732 


0.005 


9.0 x 10~ 4 


3.11 x 10 8 


16a 


18900 


250 


732 


0.006 


1.3 x 10~ 3 


3.11 x 10 8 


23a 


27200 
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TABLE VI. Average v e — * v^ and v^ — ► v^ oscillation probabilities, and muon neutrino and antineutrino CC interaction 
rates for a detector L — 7332 km downstream of a storage ring containing unpolarized muons. The event rates assume v e — •> v^ 



oscillations with sin 28 = 0.1. Also shown is the average oscillation probability for v e 



assuming no matter effects; the 



probability for u e 



without matter effects is the same as for v e 



E» 


Am 2 




NK CC) 




N(^ M CC) 


No matter effects 


(GeV) 


(eV 2 /c 4 ) 


(P(u e -> !/„)) 


(per 10 kt-yr) 


(P(P e -» ^)) 


(per 10 kt-yr) 


(J>. - !/„)) 


10 


0.002 


0.46 


90 


0.015 


1.4 


0.041 


10 


0.003 


0.53 


110 


0.019 


1.8 


0.042 


10 


0.004 


0.27 


55 


0.032 


3.0 


0.061 


10 


0.005 


0.10 


20 


0.031 


2.9 


0.047 


10 


0.006 


0.087 


18 


0.030 


2.9 


0.046 


20 


0.002 


0.11 


170 


0.011 


8.8 


0.082 


20 


0.003 


0.29 


460 


0.016 


13 


0.069 


20 


0.004 


0.46 


730 


0.015 


12 


0.041 


20 


0.005 


0.55 


870 


0.014 


11 


0.032 


20 


0.006 


0.53 


840 


0.019 


15 


0.042 


50 


0.002 


0.0067 


170 


0.0016 


21 


0.033 


50 


0.003 


0.025 


630 


0.0044 


57 


0.057 


50 


0.004 


0.059 


1480 


0.0079 


110 


0.074 


50 


0.005 


0.11 


2750 


0.011 


140 


0.082 


50 


0.006 


0.18 


4380 


0.014 


180 


0.082 
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TABLE VII. Electron neutrino CC interaction rates calculated for a distance L downstream of a storage ring containing 
unpolarized positive muons assuming v e — + v T oscillations with sin 2 26 = 0.1. Also shown are the tau neutrino CC interaction 
rates due to oscillations. 



E» 


L 


Am 2 


(P) 


N(^ e CC) 


Significance 


N(u T CC) 


(GeV) 


(km) 


(eV 2 /c 4 ) 




(per 10 kt-yr) 




(per kt-yr) 


10 


7332 








200 








10 


7332 


0.002 


0.46 


110 


8.9cr 


1.1 


10 


7332 


0.003 


0.53 


95 


11a 


1.6 


10 


7332 


0.004 


0.27 


150 


4.5cr 


0.9 


10 


7332 


0.005 


0.10 


180 


1.5cr 


0.3 


10 


7332 


0.006 


0.087 


180 


1.3a 


0.2 


20 


7332 








1580 








20 


7332 


0.002 


0.11 


1410 


4.5a 


7.2 


20 


7332 


0.003 


0.29 


1120 


14a 


14 


20 


7332 


0.004 


0.46 


850 


25a 


17 


20 


7332 


0.005 


0.55 


710 


33a 


16 


20 


7332 


0.006 


0.53 


740 


31a 


16 


50 


7332 








24100 








50 


7332 


0.002 


0.007 


24000 


1.0a 


4.2 


50 


7332 


0.003 


0.026 


23000 


3.9a 


19 


50 


7332 


0.004 


0.062 


22100 


9.8a 


51 


50 


7332 


0.005 


0.12 


20900 


19a 


100 


50 


7332 


0.006 


0.18 


19300 


31a 


170 


10 


732 








19400 








10 


732 


0.002 


0.010 


19200 


1.4a 


1.7 


10 


732 


0.003 


0.021 


19000 


3.0a 


3.7 


10 


732 


0.004 


0.033 


18800 


4.7a 


6.3 


10 


732 


0.005 


0.045 


18500 


6.4a 


9.2 


10 


732 


0.006 


0.057 


18300 


8.2a 


12 


20 


732 








156000 








20 


732 


0.002 


0.0028 


156000 


1.1a 


8.8 


20 


732 


0.003 


0.0060 


155000 


2.4a 


20 


20 


732 


0.004 


0.010 


154000 


4.0a 


34 


20 


732 


0.005 


0.015 


154000 


6.0a 


52 


20 


732 


0.006 


0.021 


153000 


8.3a 


73 


50 


732 








2.37 x 10 B 








50 


732 


0.002 


5 x 10~ 4 


2.37 x 10 6 


0.7a 


42 


50 


732 


0.003 


1 x 10~ 3 


2.37 x 10 6 


1.6a 


95 


50 


732 


0.004 


2 x 10~ 3 


2.37 x 10 6 


2.8a 


170 


50 


732 


0.005 


3 x 10~ 3 


2.36 x 10 6 


4.4a 


260 


50 


732 


0.006 


4 x 10~ 3 


2.36 x 10 6 


6.3a 


380 


250 


732 








2.73 x 10 8 








250 


732 


0.002 


2 x 10 _s 


2.73 x 10 8 


0.3a 


360 


250 


732 


0.003 


5 x 10~ 5 


2.73 x 10 8 


0.8a 


810 


250 


732 


0.004 


8 x 10~ 5 


2.73 x 10 8 


1.3a 


1400 


250 


732 


0.005 


1 x 10~ 4 


2.73 x 10 8 


2.1a 


2300 


250 


732 


0.006 


2 x 10~ 4 


2.73 x 10 8 


3.0a 


3200 
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FIG. 1. Total neutrino flux of a given flavor versus angle at the source from the beamline at L — 732 km for _E M = 10 GeV 
(solid line), 20 GeV (dashed), 50 GeV (dotted), and 250 GeV (dot-dashed), assuming 1.6 x 10 20 neutrinos in the beam per 
operational year. A Gaussian muon beam divergence of 1 mr has been folded in. The total antineutrino flux is the same. 
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FIG. 2. Total neutrino flux of a given flavor versus perpendicular distance from the beam axis at (a) L = 0.3 km, (b) 732 km, 
and (c) 7332 km (dotted) for for S M = 10 GeV (solid line), 20 GeV (dashed), 50 GeV (dotted), and 250 GeV (dot-dashed), 
assuming 1.6 x 10 20 neutrinos in the beam per operational year. A Gaussian muon beam divergence of 1 mr has been folded 
in. The total antineutrino flux is the same. 
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FIG. 3. Neutrino and antineutrino spectra for muon-type neutrinos (top) and electron-type neutrinos (bottom) in the beam 
downstream of a muon storage ring neutrino source containing unpolarized muons. 
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FIG. 4. Event energy distributions for the CC interactions of muon-type neutrinos (top) and electron-type neutrinos (bottom) 
downstream of a muon storage ring neutrino source containing unpolarized muons. 
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FIG. 5. Lepton energy spectra for CC v^ (top left), j/ M (top right), v e (bottom left), and u e (bottom right) interactions. 
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FIG. 6. The Am 2 in Eq. 17 that yields a maximum vacuum oscillation probability, shown versus the neutrino energy for 



three baseline lengths: (i) Fermilab 
Kamioka Mine (dotted line). 



Soudan (solid line), (ii) Fermilab — * Gran Sasso (broken line), and (iii) Fermilab 
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FIG. 7. Neutrino and antineutrino charged current interaction spectra for muon-type neutrinos (top) and electron-type 
neutrinos (bottom) downstream of a muon storage ring neutrino source containing unpolarized muons. The differential dis- 
tributions are shown as a function of the parameter rj = Am 2 L/E^ for rj — 0.5 (broken curves), rj = 1 (dotted curves), and 
rj = 1.5 (dot-dashed curves). The solid curves show the spectra in the absence of oscillations. 
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FIG. 8. Regions of Am 2 jE v versus sin 2 29 which have maximal mixing in matter in the Earth. The range of Am 2 /E„ values 
is due to the variation of the density in the mantle and the core. 
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FIG. 9. Regions of E v versus sin 2 29 which have maximal mixing in matter in the Earth for Am 2 = 3.5 x 10 3 eV 2 /c 4 . The 
range of E v values is due to the variation of the density in the mantle and the core. 
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FIG. 10. Density profiles along a selection of chords of length L passing through the Earth; the horizontal axis is the fraction 
of the total path length. 
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FIG. 11. Electron neutrino disappearance probability P{y e —* v x ) for x — fi or r, shown as a function of the assumed 
matter density for 10 GeV electron neutrinos propagating 7332 km through the Earth. The curves correspond to the oscillation 
parameters sin 2 29 = 1 and Am 2 as indicated. 
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FIG. 12. Modified ^ M CC interaction spectra for a 10 GeV muon storage ring neutrino source located at FNAL and a detector 
at the Gran Sasso underground laboratory, shown for several values of the oscillation parameter Am 2 , assuming sin 2 28 — 1. 
In each of the four panels the upper curves show the unmodulated spectrum, and the lower curves the modulated spectrum 
corresponding to the indicated Am 2 . 
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FIG. 13. Predicted ^v appearance CC interaction spectra for a 10 GeV muon storage ring neutrino source located at FNAL 
and a detector at the Soudan mine, shown for several values of the oscillation parameter Am 2 , assuming sin 2 29 — 1. The 
curves correspond to Am 2 = 0.02 eV 2 /c 4 (dotted), 0.015 eV 2 /c 4 (solid), 0.01 eV 2 /c 4 (broken), and 0.005 eV 2 /c 4 (dot-dashed). 
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FIG. 14. Single event per kt-yr contours calculated for v^ — > v T appearance 732 km downstream of a muon storage ring 
neutrino source. From left to right, the contours correspond to unpolarized muons with energies _E M = 250, 50, 20, and 10 GeV. 
The shaded areas correspond to the Kamiokande and Super-Kamiokande allowed region of parameter space. Also shown are 
the expected regions of sensitivity for the MINOS and K2K experiments (as labelled). 
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FIG. 15. Contours of constant v T CC interaction event rate for v^ — > v T appearance 732 km downstream of a 50 GeV muon 
storage ring neutrino source. From left to right, the contours correspond to event rates of 1, 10, 100, and 1000 per kt-yr. The 
shaded areas correspond to the Kamiokande and Super-Kamiokande allowed region of parameter space. Also shown are the 
expected regions of sensitivity for the MINOS and K2K experiments (as labelled). 
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FIG. 16. Single event per 10 kt-yr contours calculated for v E — + ^ M appearance 732 km downstream of a muon storage ring 
neutrino source. From left to right, the contours correspond to unpolarized muons with energies _B M = 250, 50, 20, and 10 GeV. 
The regions favored by the LSND results and the MSW solar neutrino oscillation solutions are also shown. 



:-!;-! 




<L.-4 



10 



-6 



MSW 



L I 



J±LLiJ I I I I U 



10~ 6 10" 5 10" 4 10" 3 10" 2 10' 
Sin 2 2i* 



FIG. 17. Contours of constant v^ CC interaction event rate for j/ e — > i/^ appearance 732 km downstream of a 250 GeV muon 
storage ring neutrino source. From left to right, the contours correspond to event rates of 1, 10, 10 2 , 10 3 , 10 4 , 10 5 and 10 6 per 
10 kt-yr. The regions favored by the LSND results and the MSW solar neutrino oscillation solutions are also shown. 
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FIG. 18. Single event per 10 kt-yr contours calculated for u e — ► i/ M appearance 7332 km downstream of a muon storage ring 
neutrino source, assuming Am 2 > 0. From left to right, the contours correspond to unpolarized muons with energies E^ — 250, 
50, 20, and 10 GeV. 
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FIG. 19. (a) Single event per 10 kt-yr contours calculated for v e — ► z^ M appearance 7332 km downstream of a muon storage 
ring neutrino source, assuming Am 2 > 0. From left to right, the contours correspond to unpolarized muons with energies 
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FIG. 20. Contours of constant v T CC interaction event rate for v e —* v T appearance 732 km downstream of a 50 GeV muon 
storage ring neutrino source. From left to right, the contours correspond to event rates of 0.1, 1, 10, and 100 per kt-yr. The 
shaded areas are excluded by Bugey and Chooz v R disappearance null results; the horizontal band indicates the range of Am 2 
suggested by atmospheric neutrino oscillations. 
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FIG. 21. Single event per kt-yr contours calculated for v e —* v T appearance 732 km downstream of a muon storage ring 
neutrino source. From left to right, the contours correspond to unpolarized muons with energies _E M = 250, 50, 20, and 10 GeV. 
The shaded areas are excluded by Bugey and Chooz Ve. disappearance null results; the horizontal band indicates the range of 
Am 2 suggested by atmospheric neutrino oscillations. 
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